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IHTRCDtJCTION 

The subcrltical assembly has been used as an educational 
device to denonstrate the principles of nuclear science and 
engineering, and for research in cxperiaental reactor 
engineering. There are several advantages in the use of a 
subcrltical assembly. The assembly is small in size and 
relatively inexpensive in aatorials, construction, and re- 
quired apparatus. Only ncminal shielding is required due to 
the low activity involved. Lattice parameters can be varied 
more easily than in a critical assembly. 

The kinetic equations of the assembly are tine dependent, 
non-linear differential equations which nay be linearized in 
the derivation of the transfer function. The method of ex- 
perimental evaluation of the transfer function is to vary one 
or nor© of the parameters in such a way that the limits within 
which the linearisation is valid are not exceeded. This 
variation of the parameters is called "forcing”. The usual 
methods of forcing in a subcrltical assembly are varying the 
strength of the neutron source or the reactivity in a periodic 
mariner. The minimum periodic input signal to which the re- 
sponse can be ncajjured is determined by the amplitude of the 
random or stochastic variations in the steady state flux 
level due to the inherent statistical nature of the fission 
process. It is possible, in a system with large stochastic 
variations in the neutron flux, to exceed the range of line- 
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arity of the kinetic equations fey making the magnitude of the 
input signal largo enough to obtain a readable response. A 
technique has been developed, from studies of stochastic 
processes such as Brownian Motion, which makes it possible 
to evaluate the transfer function of an assembly. The 
stochastic or random variation in the flux is considered to 
bo the response of tho system to stochastic source forcing 
duo to the statistical nature of the basic phenomena involved. 
This method of evaluation of the transfer function should be 
least likely to invalidate the linear approximations made. 
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EEVmi 0? THF LITFEi\TUBF 

A great anoimt of effort has been expended in the deriva- 
tion and evaluotlon of critical assembly transfer fnnetions. 
One of the primary reasons was the need for an autcastic 
control system to mintain a given pov;er level in a reactor. 
The similarity of the transfer function of a reactor to that 
of a servomechanism control system had been noted, and it was 
determined that this type of system was suitable for reactor 
control. 

Ov/ens, Crever, and Pigott (13 ) in 19^9 proposed the con- 
cept of treating a critical assembly as a "black box" to 
which a feedback control system could be applied in order to 
maintain a desired power level. The proposal was based on 
the observation that the roots of the kinetic equation of the 
assembly Moro negative. This indicated that the critical 
assembly transfer function was the same as the transfer func- 
tion for a minimum phase shift network, and that the transfer 
characteristic could be represented by a Bode diagram. In 
that same year Frans (!?) presented a more detailed derivation 
of the transfer function and showed that the standard tech- 
niques of servomechanism theory applied. Harror, Boyar, and 
Krucoff (9) evaluated the transfer function of the CP-2 re- 
actor by sinusoidally varying the reactivity, and showed that 
the transfer function could be used in the design of servo 
loops or regulating systems which include the reactor. 
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Mockup assemblies have been utilized for the determina- 
tion of the transfer function of proposed reactors. Boland j 
Smith, and Kice (3) measured the zero po\^er transfer function 
of a fast critical assembly, 2PR-III nockup assembly of 
EBR-I, Hark III, b/ determining the change in power level due 
to oscillation of the reactivity. 

heflnemont in the techniques of testing and the evalua- 
tion of transfer functions are constantly sought, and the 
subcritical assembly provides a relatively simple vehicle 
for such Investigations. A sxibcritical assembly was utilized 
by Aztraann, Dessauer, and Parkinson (1) for reactivity meas- 
urements and testing of pile components. The transfer func- 
tion of the Iowa State University subcritical assembly Mas 
Initially evaluated by Fdcci (l4) utilizing the method of 
source forcing. 

A new technique for evaluation of the transfer function 
of a nuclear reactor operating at a steady state by an analyst 
of the random, or stochastic, nature of the neutron population 
has been developed recently. Koore (11) has shewn mathemat- 
ically that the square of the modulus of the transfer function 
of a reactor is proportional to the Fourier transform of the 
autocorrelation function which is equal to the power spectral 
density function of the system. This was verified experi- 
mentally by Cohn (4) and by Griffin and Lundholia (8). Velez 
(16) evaluated the transfer function of the Ford Nuclear 
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Reactor and tho aiibcriticval ^sscncly at the University of 
ilichigan by the direct measurement of the autocorrelation 
function, but the results were reported to be inconclusive in 
the case of the subcrltical assembly. All references found 
describing the evaluation of the transfer function by measure 
ment of the power spectral density function pertained to 
critical assemblies. 
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'Ihe kinetic equations of the siibcritlcal assembly can 
be derived from the general diffusion equation (6, p, 223 
7, p. 101) 



v;hich is applicable to that region of a subcritical assembly 
away from the extraneous neutron source where most of the 
neutrons are fission neutrons. The syrabcls used are defined 
in the table of symbols. The first term represents tho less 
of thermal neutrons due to leakage from the assembly; tho 
second tern is the loss of thermal neutrons due to absorption; 
the third term is the neutron production or source term; and 
the right side cf tho equation is the non- equilibrium changc 
in neutron population '-’ith time. 

I'he source term is made up of prompt and delayed fission 
neutrons, and neutrons from an external source. It has been 
shown (7, p. 193) that the number of fast fission neutrens 
produced is 0, and further (6, p. 226) that if ^ 

is the total fraction of delayed neutrons the prompt neutron 
source term is 



- Z^0 + S 

iX 



Zn _ 1 ^0 

^t V ^ t 



( 1 ) 




( 2 ) 



There are six laiown groups of delayed neutrons that 
account for 0.64 per cent cf the thermal neutron production 



V 



M H •••••vmip mt tawiCMM •UM«f «lt 

wm.m --r ' '!-» j 

f J^Ult ta «i4**CA|f|t fJ Riflll 

■lif It f4i» 4««t% «TM RTfM €MMM tmn ymf 

m tiRprif i i MM 9tm <iir .aJiant %• •« ar 

ttt nUMM ib RC1 Miii tt R« b 

«# i^ tRiRa t X4R»4^ U tUI mii H «*««i mrr 

iMt \mm mmm ^ «i 4MM 

at>ttA2Jtt*R« tv ti RAiMt tV tt #l4t *4^H ^ 

^ m>. 

mWitf VK Hflii t At t Ml attlMb 
cmHaH VRi£»i 4wi#p«t| ft ^ ^ ttVAt Mtt 

t»t ^1 .«*»•• fMtfa* m ^ Rfft it R .tR ,mmh 

*\u «•« Itee VI 4« ,%~\ 

»9m0% mU i ii ti i i tt !• tWR? Uitf tV U 

•I 1^ VM 




7 



These delayed neutrons accompany the radioactive decay of cer- 
tain fission products, called neutron precursors, which exist 
in an excited state. the fraction of each delayed 

neutron group then the formation of the precursors of the 
group would be V.’hen the concentration of a precursor 

group is denoted by Cj|^ and the decay constant the rate 

of decay would be/^^c^. The net rate of formation of pre- 
cursors to the 1-th delayed neutron group would then be 




The total rate of production of delayed neutrons is ^ Xc.. « 

i=l ^ ^ 

and when multiplied by the resonance escape probability, p, 

-B^ T 

and the slowing down non leakage probability, e , the 
result is the delayed thermal neutron source term 

S. = p ZI \ c. (U) 

° i=l ^ ^ 



Table 1. Table of symbols 



Buckling 

^ Fraction of delayed neutrons 
^ ^ Fraction of 1-th group of delayed neutrons 
C Number of delayed neutron precursors in the 
assembly 

Number of the 1-th group of delayed neutron 
precursors in the assembly 
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Table 1. (Contirue')) 



c 

D 

f 

k 

k<» 

L 

r 

io 

I 

A 

H 

n 

p 

0 



S 



e 



Concentration of delayed neutron precursors 
Concentration of i-th group of delayed 
neutron precursors 

Thermal neutron diffusion coefficient 
Frequency 

Effective multiplication factor 
Multiplication factor for an infinite assembly 
Neutron diffusion length in the moderator 
Prompt neutron lifetime in a finite assembly 
Prompt neutron lifetime in an infinite assembly 
Average neutron lifetime in a finite assembly 
Average decay constant of delayed neutron 
precursors 

Decay constant of the i-th group of delayed 
neutron precursors 

Total number of theriaal neutrons in a finite 
assembly 

Number of thermal neutrons per unit volume 
Resonance escape probability 
Therraal neutron flux 

fate of production of therml neutrons per 
unit vclume from all source 
Total external th -rnal neutron source 
contribution to the assembly 



cm. 

cycles 

sec. 



cm. 

sec. 

sec. 

sec. 

see. 

sec. 



-1 



-1 



neutrons 

3 



CIS.' 



neutrons 
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cm. sec. 

neutrons 
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cni.“^sec. 



neutrons 

sec. 
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Table 1. (Continued) 



External thermal neutron source contribution 





per unit volume 


neutrons 

3 

cm.-^sec. 


2^3 


Macroscopic absorption cross section 


cm. 


t 


Time 


sec. 


nr 


Fermi age 


2 

cm. 


V 


Thermal neutron mean velocity 


cm. /sec. 


CjO 


Frequency 


radians 
sec . 



The contribution of the external thermal neutron source term 
will be assumed to be 



S = ^ 



(5) 



Substitution of Equations 2, 4, and 5 into Equation 1 yields 



■|£ = D^20 . + (1 )k, z^0e 

+ p e"® ^ZT L 0, + 

i=l ^ ^ 



-YT T 



( 6 ) 



0 and c^ are functions of space and time, and it can be shewn 
(6, p. 227) that the space and time variables are separable 
allowing the partial differential equation to be reduced to 
an ordinary differential equation. Consideration of the 
fundamental space mode only (6, p, 227 and 9, p. 33) yields 
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the eouallty 



V ^-0 



= - B ‘^0 



(7) 



The diffusion length is defined by 






( 8 ) 






Eearrangenent of the equation yields 



D = L^Z 



'B. 



(9) 



The aean life of thoiTial neutrons In an infinite nediun. is 
defined by 

_ 1 



O = V " ZgV 

When rearranged this relation becoraos 



(10) 



z =-^ 

a / V 

-^o 



( 11 ) 



The flux, the effcctivo ml ti plication factor in a finite 
aediu!?!, and the pronpt neutron lifetine in a finite mediun 
are defined by 

0 = nv (12) 
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(lb) 



The combination of equations 13 and lb yields 
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A k 



( 15 ) 



The substitution of I-quations 7, 9j 11 j 12, and 15 into I'qua* 
tion 6 yields 



H = - + 1) ^ + (1 



( 16 ) 









The combination of Equation l4 with iquation 16 and rearrang- 
ing yields 



(17) 



2 ° 

^ = [k(l -/^*)- l] ^ 'I- p ^Z A.c>+ 
dt / / i--=l 

When Equations 11, 12, and 15 si's substituted into ijuation 

3 the result is 



dc 



dt 



1= /t 



k n 






" ^i^i 



(18) 



pe 



Equations 17 and 18 are the coupled Vrinetic equations of the 
sufccritical assembly. 

It should be noted at this point that in arriving at 
Equation 17 the assumption has been made that the external 
source contributes a different number of thermal monoenergetic 
neutrons to each unit volume of the subcritical assembly, 
the number decreasing with height in the assembly. The 
kinetic equations can also be derived by considering the 
total neutron population rather than the neutron population 
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on a per unit volirae basis. 

Consider that at a particular time the assembly contains 
N neiitrcns, and the effective reprcducticn constant, k, is 
less than one. 'Ihe prompt neutron lifetime is J/* ^ the prompt 
neutron reproduction constant is k^, and the rate of decrease 
of neutron population in the assembly is (k^ - 1) if 

neutrons are not supplied to the assembly from an external 
source. Absorption ana leakage account for the decrease. 

If is the fraction of delayed neutrons then the prompt 
neutron re-production constant, k^, can be replaced by 
k(l -y^). The rate of decrease of neutron population uith 
time is then ~lven by 

If = [k(l -/? ) - l] U// (19) 

Neutrons are supplied to the assembly by an external source 

and six groups of delayed neutrons. The total number of 

delayed neutrons of the i-th group stored in the fissions 

products of the assembly at any time is If the decay 

constant of the nuclei containing the i-th group of delayed 

neutrons is the neutrens of that group are emitted at 

a rate as fast neutrons. These fast delayed neutrons 

T 

must be corrected for fast r.on-leakage, o” , and resonance 
escape probability, p. Sutnmaticn of the terms in the non- 
equilibrium rate equation for the assembly yields 
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dt 



[k(l l] ^ (20) 



The total number of thcrnal neutrons being produced per 
unit tiao is of which a fraction are delayed. Ihe 
delayed neutrons are not lost by leakage or absorption until 
after they are emitted from the fission products, and hence 
the number of delayed neutrons stored in the fission products 
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is greater than by the factor 
rlum equation for the i-th group of delayed neutrons is 



d C 
”dt 









i '=1 



< 21 ) 



Equations 20 and 21 are the sane as Equations 17 and 18 ex- 
cept that they aro based upon the total number of neutrons 
in the assembly, and the external source term then accounts 
for all of the neutrons supplied to the assembly. The neutron 
flux at any point in the assembly can be considered to be 
proportional to the total number of neutrons in the assembly 
at a X)articular instant. The use of cither set of coupled 
kinetic equations in the derivation of the transfer function 
will yield the same result. 







'WM ^ n«r 



fMMUfMM»n wa^M maaH mmm^^ 

Vmmimms^Mki y nim il# fgl •^•nmrS ^u 



fl»|i 





f 




i» mma •ds #m 93^ ti*» ^ Mii«i^ 
«Aw Ml «rf* 4 ms^ rv« #fedf Ipw 

<0«s MlMM'liasvAtM ^^4 b# Hf 6i 

iw«« <1 ^iiv^ wfm»9m «#i !i« HH 

4 Ml #< «nrTfif»fM «# di amar^ ^ mt% 



la 




i *vMT# %#M4irr' 
•# n 9 % ta lTi ^ W ii 




hua% liibi 



df I 



DEBIVATION OF IKI TRANSPH.. FUMCTICHS 



Phe kinetic €qu'’.ticns and the reanltant transfer function 
can be varied by varying the scurce or reactivity in a pre- 
detoz”^ ined manner. This variation is usualJ^y sinusoidal when 
purposely intrcdluced. Howci/cr, the scurce strength and re- 
activity raay vary even at "steady state" conditions. This 
is due to the fluctuation in the number of neutrons available 
to the assembly caused by the natural statistical fluctua- 
tions in the rates of neutron absorption and fission. The 
resultant randem variations in the neutron flux make possible 
the evaluation of the transfer function by proper analysis 
at steady state. 

Transfer Function with Source Variation 



For steady state cenditions (using zero subscripts) 
Lquations 17 and 18 become 

6 

( 22 ) 



an 

dt 



= 0 = [k(l -^)- l] ^ pe 






^ A . c . + 
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-b^t 7 
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i ^i. 



(23) 



Consider a snail variation to be superimposed on 
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sc that ^ ^ c 

n = ^ n 



(24) 
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Substitution of Lquation 2^ into Fquations 17 and 18 yields^ 



dt 



A r -1 n) 

^ = [k(i -/ ) - 1] — 

6 ^ 



dt 



p© 
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and 



dcj d < 



dt dt 



j /^k (n. + <^ n) >. ^ 

i ^ —P — ) (c^ ^ cTe^) 

pe“^ ^ r i i 



(25) 



(26) 



When the tor/ss of fquations 25 and 26 arc expanded, the 
steady state compononts of i'quations 22 and 23 arc collected 
and set equal tc sero, the resultant equations are 

^ n pc“® ^ ^Z X^ S c^+ (27) 



d cT ( 



dt 



. -d 2^ n - /jCj 






(se) 



hy transforiaing to Laplace notation Equations 27 and 28 be-coniG 

2 ^ 

s^n(s) a c^n(s) - ^n(s) + pe“® ^ ZZ \.fc^(s) 

A A i=:l 



+ ( s ) 

s<^c^(3) = <rn(s) - Xj <^c^(s) 

pI e 

Rearrangement of Equation 30 yields 

r , ^ fyi 1 r „ 



(29) 

(30) 
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Substitution of Equation 31 in Equation 29 and collecting 
terms yields 

6 

cTnCs) [ s - ^ J= = > (32) 

By rearranging Iquation 32 ^ the transfer function becomes 



c^n 
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(33) 



SJ’ - k -M -M5^ - k 



By substitution of 



S + A 



^ A h 

i’ot’ Z_ trr — ^ Equation 33 becomes 
^ i=l ^ '^i 



cT 



J' 



s/- k(l -/’) + 1 - k 



(3^) 



Hultlj)li cation of the numerator and denominator by (s +\) 



and division by J yields 



44 = 

^(s) 



A 



S2 + ( X + ^)g + ^ J,j 



(3?) 



The denominator is of quadratic fcrni, and when average values 

of A = 0.08, k = 0.5, J* ~ 0.001, and = 0.006U are 

substituted in the coefficient of s, it is seen that A and 
^ 1 - k 



are negligible compared to — so that 
cAn(s ) s •<• A 

- (^)s + A-ii^ 



(36) 



The roots of the denominator are found by use of the quadratic 
formula . 
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= ^ -li 



' t*Ai* 
1 - ITT 



From the binomial expansion (a + b)^ = a where 

r U A o\ 

b 3) the radical expands to ll - - 1 - ^" " 1 "' ^ 

By substitution of the expanded terms for the radical, s 
becomes 
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1 - k 



First root 

s - l_=Jc 

®1 2 J* 

Second root 






s - k L 

5o - o /; ♦ T - 



, - 

’2 2/' 



1 - 1 + 
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] = - + A 



1 - k J - ' 

When A is considered negligible comparcxi to the first term 



the resultant root is 
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The denominator is now seen to bo factorable to a good ap- 
proximation into 

(3 + A )(s + 

By substitution back into Equation 36 the result is 



n(s) 



j,, t,„A 1 

(s + X )(s + j¥^) s 



( 37 ) 



This is the final form of the transfer function. It shews 
that for a subcritical assembly one may ignore the effect of 
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the delayed neutrons in determining the change in neutron 
population due to a variation in the source. The terms con- 
taining the average delayed neutron decay constant cancel 
out, and thus the prompt neutron lifetime can be replaced by 
the average neutron lifetime. The assumption that the effect 
of the delayed neutrons can be ignored is valid for 
0,995 (l)* ihe subcritical assembly utilized in this in- 
vestigation has a ot about 0.52, hence the above assump- 

tion should be valid. The same result, i-quation 37> is ob- 
tained from Iquations 1? and l8 when the delayed neutron terms 
containing are deleted and the same derivation procedure 

as above followed. 

Transfer Fxmetien with hoactlvity Variation 

Again utilizing Equation 17 and Equation 18 for the 
transfer function and '•-quatiens 22 and 23 for the steady 
state condition, let a snail variation k be imposed upon 
k so that 

k = ko + c^k 

n = n^ + (38) 

Cl = Ci^ H- <fci 

where cTk is the amplitude of the variation, and not the 
instantaneous value of the excess reactivity. Substitution 
of Equation 38 into Lquations 17 and 18 yields 
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+ pe-®^'^ E A,(c, + i-c.) 

i=l ^ 



/* 

'4' -'^ 



(39) 



dc4 d<^c. _ (k - <^k)(n + cf n) . 

- ^ 7 - ^ A,(c (40) 

p / e-® ^ ® 



-JL 

dt 



By multiplication, setting products of small numbers equal 
to zero, and collecting the steady state terms of Equations 



22 and : 


23 and 


equating them 


to zero 


one obtains 




d n 
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, "o'''" /j"o 
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n^cTk 




dt 
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n 
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2 6 

+ pe 

i=l ^ 


c^Ci 




(41) 


ifii 

dt 


^2 T- ♦ ra2 rp 1 

pi’e'*’ ^ p(>e-® ^ 


cT 


(42) 



By transformation to Laplace notation the equations become 
s<Tn(s) = <^k(s) + ~ cTn(s) - cTn(s) ~ ° cTk(s) 

- ^ + pe”^ ^ZI cTc. (s) (43) 

i* 1=1 ^ ^ 

(3) = '^n(s) + cTk(s) “A4 cTc4(3) (44) 

p/e p/e 

Solution of Equation 4U for cTc^(s), substitution in Equation 
43, and rearranging yiolds the transfer function 
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The substitution of for ,ZI ^ in Fquation 4-5 

1 **i i 

yields 
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S + A 

The multiplication of numerator and dcncminatcr by s + /^ and 
collection of terms yields 
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Division of the numerator and denominator by J*/k^ yields 
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(48) 
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This equation is seen to have the same denominator as Iquation 

1-k 

35 which simplifies and factors to (s +A)(s + ) . In the 

numerator the factor (1 -/?) is seen to be nearly equal to 
one when 0.0064 is substituted Tor ^ . Fcuation 48 then be- 
comes 
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which is the final form of the transfer function. This result 
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is the same as Equation 32 where the source strength was 
varied except for the constant multiplier n^/ The a\'erage 
decay constant A and fraction of delayed neutrons drop out 
of the final expression as before, shcv;lng that the transfer 
function for a subcrltical assembly is independent of the 
delayed neutrons and dependent on the average neutron life- 
time. 



Modulus of the Transfer Function 



1-k, 



Division of the transfer function, Iquation *+9, by — -jp" 
and substitution of 3 = , yields 
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The rest It of noraaliaing Iquation 50 and setting it equal 
to co) Is 
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It is seen from Ii:iuotion 5l that the prompt neutron break 
frequency, co^, is equal to for a subcrltical assembly. 
6Dp eff€-!ctively determines k for the subcrltical assembly if 
J is known, and the reverse is also true. The numerator 
and denominator sre multiplied by the complex conjugate, 
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The modulus of the transfer function, or the nagnitude of 
Yi^co), is 







which reduces to 



|y(j«)| = 






(53) 
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The square of the modulus of the transfer function 
is then 
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The square of the modulus of the transfer function is the 
basis upon which the parameters of the transfer function can 
bo evaluated as will be shewn in the next section. 



** ! 



f 







It r. 






§ i » i ^ 



f - 

i 











\J *...*« ^ 



ipy 

*' ji- r t^l gu|)ti 



't>/ 



%X 



r^' • 









c 



§n<i-4 jQa >? ifl' 

-»ilf u 

— • - •!t«i)»| 



Tifc 



•it I ••€?•<'• ^ ••If \i| Ujf? 



mmi irri^mrqf^ r*'^«i^*t# * "**f »•!< ^14# cofff s^i%g 9 

.^0j^fmek 9cmm ••# pat •• lli» •» «4 



e 






23 



STOCHASTIC rVALUATION' 0? THT TEAK3FEK FUNCTION 

The transfer function of an asseiably can be detorained 
exp crlsien tally by introducing sinusoidal variations in the 
input and measuring the phase angle and amplitude attenuation 
with frequency of the resulting sinusoids! variation in the 
output. This is usually accomplished by sinusoidally varying 
the strength of the neutron source or the reactivity and 
noting the resulting x’ariation in the neutron flux. Vari- 
ations in the neutron population due to the statistical 
nature of the fission process cause a change in reactivity 
and a fluctuation of the power level of a critical assembly. 
These random variations in the flux level nay be considered 
as a suiiiraation of Fourier series which transfora into Fourier 
integrals whisn the period approaches Infinity (2). The 
transfer function of an assembly could then logically be ex- 
tended to represent the response of a system to a random as 
well as sinusoidal input. 

Consider a voltage which is a random or stochastic func- 
tion of time and limited to a small range of frequencies be- 
tv/een f^^ and fp. On© can detenrino a value which is propor- 
tional to the average power by squaring the function, inte- 
grating, and dividing the result by the time over which the 
integration xms performed. If this "av^'rage power" is divided 
by Af - f]^- fpj the frequency increment which the average 
power represents, the result is the average power or mean 
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square voltage (amplitude) per unit frequency. This is called 
the pover spectral density function. The term power spectral 
density function has been extended to apply to any function 
which consists of mean square signal amplitude per unit 
frequency of the band width investigated. 

Nuclear reactions and fission processes which give rise 
to the ntsutron flux in a subcritical assembly can be con- 
sidered as a random variable which has a normal or Gaussian 
distribution about some mean value of frequency of occurrence. 
The power spectral density of such a random variable is the 
average power (mean square signal amplitude) per unit 
frequency of tho band \ddth measured (2). It has the unique 
property of being a constant for this typo of function and 
is defined as a "white noise". When such a random variable 
is passed through a linear system the output is a random 
variable which has been attenuated by the system character- 
istics (10). 

If a system such as the subcritical assembly is excited 
by a white ncise, which has a constant amplitude at all 
frequencies^ one can analyze the output of the assembly over 
a range of frequencies to determine the response or transfer 
function (8). It has been shown mathematically that the 
output power density spectrum of a chain reacting system is 
proportional to the square of the modulus of the transfer 
function (11). This has also been demonstrated experimentally 
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for several nuclear reactors (4, 8). 

Figure 1 shews on if?cal cysten for power spectral density 
measurement in which it is assumed that the electronic 
components have perfect frequency response over the range of 
the investigation and introduce no ncise. Under these condi- 
tions the system equation would ho 

Gyy(W) = G^( w)|Y(3i^)l2= K|y(jco)|2 (55) 

where G ( ) = input power spectral density = constant. 

G (w) = output power spectral density. 

= system transfer function, 

1^(3 = square of modulus of transfer function. 

Since the frequency response of the measuring equipment 
is net perfect, it alsc introduces a transfer function into 
the equation. This nay be neglected and the response of the 
measuring equipment considered independont of frequency, if 
the break frequency of the oquip’~ent is beyond the range being 
investigated (8). The final form of the system equation was 
shown to be (6) 

Gyy( w ) = A + b|Y(J^o )|2 (56) 

where A represents the level of measuring equipment noise 
and B represents the level of the assembly noise. The sub- 
stitution of Iquation 5^ into Lquation 56 yields 
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Figure X. 



Ideal systora for dotcrreln'itlon of output rowr 
spectral density function 
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G ( ^ ) = A + 



B 



■yy' ' " 1 + 



A sketch of Iquition 57 is shewn in Figure 2. 
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Figure 2* Idodl pcwer spectral d^.slty function of an 
egsenfcly 



DESCRIPTIOH OF APPARATUS AND rXprpmr 2 «’TAL PI^OCLDITRE 



The graphite moderated subcritical assembly (lU) utilized 
in this investigation was constructed at Iowa State University 
in 1957 . Figure 3 is a photograph of the suhcritical assembly 
with the end cover removed. The blocks were constructed from 
AGR grade graphite cylinders with diameters of 7*0 in. and 
6 3/8 in. , and length of 60 in. The lower nine layers were 
nad© up of blocks with a 6.0 inch square cross section and 
rounded corners of radius 3*5 inches. The top five layers 
were constructed of blocks with a 6.0 inch by 5.0 inch cross 
section with rounded corners of radius 3 3 /l 6 inches. The 
density of the graphite was 1.56 gn. per cu. cn. or 97*3 Ih. 
per cu. ft. An 8^ inch lattice was formed by placing uranium 
slugs in alteiTiat© channels as shown in Figure 3 . The slugs 
of uranium were 1.000 in. in diameter and 8.00 in. in length 
and encased in aluminum cylinders O.OhO in. thick and v/ith 
ends 0.20 in. thick. The outside dimensions were 1.080 in. 
in diameter and o .40 in, in length. The uranium density was 
19.0 gm. per cu. cm. or II86 lb. per cu. ft. 

The assembly v/as covered on the top and sides v^lth a 
0.010 in. sheet of cadmium sandwiched between a 3/8 in. sheet 
of plyv/ood and a 1/8 in. sheet of masonite. A11 effective 
"black wall" for thermal neutrons was produced by the cadmium. 
The horizontal distances between the side sheets of cadmium 
was 60.5 in. and 62.5 in. With the extrapolation lengths 



«t » — i n n «Mf at pm* •pm «r 

^ mJ i«f ip c JiA •vtlclln «Uar«rf 4M 

s*^i «ci 14 litlftl Ittp ^•«| tV|*4 

^4i^» vvcYv* mimtm iMf 4ii) • 4ll» fV^Xi ^ •4M 

4ii^C4n» #4# W9 •4k«a4il |«| mtHm !• 

^mprn m m^l ^*4 S ill« y run 

•¥^ «Mua ilM f 44IMI % •MM* MMl tflta j»4t%nb 
4%t f4 <4^ «•« «4€ ••% if •#» UUivm •‘IK^ 1» Hlvyi 
giWft IPtnld # ftMC 4«af0/ 4M1 |f 0& lit dii -Xf 
wJ# «l amtf 4^ mmi ii •<■■1^ «4 wSi 

4hB4X aM .#1 4«4 ffKiWiI ft ^ M9,| nt^ MttMf » 

« 

mtM ^ rftAci/11 Miiinu mi ^%mm ^ 
•41 c4«i: Mmm 4wi4Mr^ tMnt|%A »4 d 14* #ia ^#9 piv* 
••f m^mm «fi ^mma ti «4f #*rf w u$pppU m 

itfi ^ m m ^ ^ ^ ^*t4 

a «■•» M|« #41 %4i IP tPi fiPmjmM pMt 

T«yb *«FJ a pp/mHMpm pm^pppa ^ tam «a4 41^4 

«a i4»Pip«i itf #«44i vit 4 iPt JUCMf^ It 
■ i^lfi^M mam fl f aa wn apa toittax CPM4r i4l * XJpt 

mymm £t mPpp^p pM naaalaf mmP^miA J«lc«cHfl# tti^ 
o#t^ i ^ifaSp^arOm^ pdl 4114 «t| ftH 4 p» miI f 44 4it 



30 






«i0 ^ • 






i 



Figure 3* Front view of the subcritical assenhly 
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added, the effective nlzo of the assenbly was approxinately 
62 in. (side to side) by 64 in. (front to back). 

The subcritica?s. assembly was mounted on a pedestal as 
shot.!! in Figure 3* ^he throe spaces under the assembly housed 
11.5 in. deep water tanks. The water provided protection 
fren the neutron scurce which was mounted in the center tank 
under the pile. 

The neutron scurce consisted of five individual 

plutcniuffi-berylliun sources in containers 1 in. in diameter 

and 1 3/8 in. in length. The individual scurce strength was 

rated at one curie emitting approximately I .63 x 10° neutrons 

6 

per sec. The total scurce strength was 8.15 x 10 neutrons 
per sec. The five sources were mounted on a pedestal in the 
center tank which was flll»='d with water. The sources v/ere 
located on the centerline of the assembly directly belov/ the 
bottom layer of graphite. Five holes were cut along the 
center line of the front panel of the assembly at one foot 
vertical intervals. The holes were cut to correspond to 
vacant channels in the lattice for insertion of the BF^ 
neutron detector. The lox^est hole, one foot above the bottom 
of the assembly, will be referred to as hole one, the hole 
two feet in height as hole two and so on. 

The electronic circuitry for meastiring the random varia- 
tions in the neutron flux is shown in Figure 4. Tho neutron 
probe used v/as a Wood Counter Laboratory BF^ detector. A 
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Fipn*© U, Eicctponie circwltry 



3 ^^ 

Nuclear Chicaco ITLtraacalsr Model 192A was used to supply 
high voltage tc the dotectcr and to aiaplify the signal fro?a 
the detector, k Tektronix Type 112 Direct Coupled Amplifier 
uith a flat frequency response from DC to 2 megacycles v/as 
used to further amplify the signal which was fed into a 
Khron-Hite ultra lov frequency Band-Pass Filter, Model 330 A. 
The band width of the filter was adjustable between a low 
cut off frequency of 0.02 cps tc a high cut off frequency of 
2000 cps. Frequency adjustment was accomplished by means of 
two dials calibrated with a logarithmic scale reading direct- 
ly in cps from 2 to 20 and multiplying switches of 1/100, 

1/10, 1, 10 and 100. A capacitance filter in the input 
circuit filtered out any DC component of the signal. The 
filter comprised an extremely low noise circuit, with In- 
ternally generated hum and ncisx? less than 100 microvolts. 

The filtered rand cm signal v;as then fed into a two channel 
Brush Kecordcr Hark II for further amplification and record- 
ing. The recording pen sensitivity was variable from C.Ol 
to 10 volts per chart line, and the chart speed was variable 
from one tc 125 ma per sec. This circuitry was used to in- 
vestigate the frequency range of 0.02 cps to 20 cps. An 
alternate system v/as used for the frequency range of 20 cps 
to 100 cps and beyond, to facilitate data taking and 
processing. All of the components after the Nuclear Chicago 
Ultrascalar v;ere replaced with a Hewlett-Packard Wave analyser 
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Model Th'S •.mv-' nnalyzer h-id an adjustable band pass 

of 30 cps to 1 U 5 cp3> and a frequency range adjustable froa 
aero cps tc I 6 kc. Ihc output ?-ieter was graduated in V'Ao 
nllllvolts. Figure 5 is a photograph of the equipment 
described above. 

Data were taken in each of the five holes, with the 
probe on the centerline of the assembly. The frequency range 
investigated was 0,02 cps to 100 cps. The circuitry including 
the band pass filtc’' and recorder v;as used for the range of 
0.02 cps to 20 cps. The wave analyser was used from 20 cps 
to 100 cps. In order to minimise the introduction of spurious 
signals into the system, the output was taken from the scalar 
without the counting circuit being activated. The signal was 
then fed into the Tektronix Type 112 anplificr xvhorc the 
input and output slgnai was chocked periodically with an 
oscilloscope to insure rdninur. a.c. pickup atid no saturation 
of the .amplifier. The signal was then passed through the 
band pass filter at a constant per cent band pass such as 
0.02 cps to 0.04 cps, 0.04 cps to O.OS cps, 0.08 cps to O.I 6 
cps, etc. The output of tho band pass filter was recorded 
on both channels of the Brush recorder for two minutes for 
each band pass width setting. Typical random data for 
representative band pass widths arc shown in Figure 6. 

The power spectral density function for the 0.02 cps to 
20 cps range was obtained for each band pass width by a method 
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Figure 6. Typical random data for representative band-pass 
widths 
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f = 0.0375 cpt 




4f«005 cpi 
? *0.07 5 cp s 




fsO.15 cpt 




4f * 0.20 c p s 
f*0.30 cpt 




All chart tpetdt 5 mm p«r tec 



id t s| 0 sec 



similar tc thot outlined by Ilurphy (12), One hundred 

ordinates of c?ach randen trace were measured at one second 

intervals. Kach ordinate was squared., and then the ordinates 

and their squares were summed separately and avenged. The 
2 

variance, 0~ ^ was obtained by taking the .difference between 
the average of the squared value and the square of the average 
value. This is equivalent to the average of the square of 
the difference from the mean. 

CT"^ = (y2 - y2) = (y- y)2 

The power spectral density function G„y( ) was then obtained 
by dividing the variance by the band pass width in cps. This 
resulted in a value of mean square amplitude per unit frequen- 
cy, or mean power per unit frequency, for the band width in- 
vestigated. 

O 

G (i^) s mean square airrolltude 

yy'- ^ Af^ unit frequency 

The power spectral density for the 20 cps to 100 cps 
range was obtained by squaring the reading from the wave 
analyzer output meter. The meter read directly in PiJS power 
per unit frequency of the band width investigated. By squar- 
ing the reading, the mean square amplitude or power per unit 
frequency of the band width investigated was obtained. 

A sample calculation of the type used to reduce the data 
from the Brush recorder is included in the Appendix. 
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DI3CU3SI0N cr rCSULTS 

The results of the data reduction are presented 
graphically in Figures 7 and 8 with the power spectral densi- 
ty, G (dj), in arbitr-’ry units versus frequency in cycles 

J ^ 

per second. The data v/sre plotted only frotn 0.03 cps to 
10.0 cps because the power spectral density appeared to bc- 
coae constant by 10.0 cps and renained a constant frera 10.0 
cps to 100 cps. This constant level represents the "A” terra 
in Equation 57 which is the white noise of the inotrucienta- 
tion which is due primarily to the DF^ detection chamber. 

The instnrnentation noise decreases vdth increasing height 
In the assembly. This is explained by the fact that most of 
the instrumentation noise is due tc the flux level or count- 
ing rate in the detection chamber which decreases with in- 
creasing height. 

Figure 7 shows the change in power spectral density 
function along the axis of the assembly. For hole one (s = 

1 ft) the power spectral density was determined tc bo ap- 
proximately a constant throughout the range of frequencies 
investigated. This indicated that the rsnden variations in 
the neutron flux were primarily due to the flux from the 
source which is "white noise". Hole two (z = 2 ft) shows a 
deviation from the constant power spectral density function 
of the instrumentation and indicates an increasing contribu- 
tion from the fission flux and external source which have 
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siiNfi AavHiiaav [(m) AAo] Noiiownd AiisNaanvaioadS aaMOd 



FREQUENCY Cf 3 - CYCLES/SECOND 

Figure 7» Comparison of power spectral density function with height in 
the subcritical assembly 



^3 



> 




Figure 8. Comparison of measured with calculated power 
spectral density function 



been attenuited by th ? transfer fnr^ticn of the asser.bly. 

Holi? three (z - 3 ft) is 91*^^ cs o* approxinatoly two diffusion 
lengths atov-: the cecternal source. 'Ihc powsr spectral density 
fraction of th : neutron flux at tf is level indicates that the 
transfer funetion of the assembly has operated on the flux 
to a nuch ;^reat r extent. Hole four (a = ft) shews a 
Etaxinuai slope or attenuation of the pov/er spectral density 
by the transfer function of the assembly. Hole five (a = 9 
ft) shows a decreasing slope and a tendency to approach a 
break frequency toward the lower frequencies. This night be 
explained by a decreisc in the ”'-f' coofficiont of Iquation 
57 vdiich represents th>- level of assc’-ibly ncisa snd is pro- 
portional tc th*‘- square root of the powf^r or flux level in 
the assciably. ks ''b" decreases, tho inotru"’'‘'ntation noise 
level app’r’oacheo th.i asscnbly ncin- level, and the rounding 
off of the breaks in th-: curve t"'nds to -aask the true slope 
of thA power spectra' d^,nsity function, 

Ko definite break frequency uas discernible frox. the 
data presented in Figure 7- The data for holes three, four, 
and five were replotted on Figure along v’ith a plot of the 
theoretical squared modulus of the transfer function from 
Iquatlon 5*i‘‘ The theoretical equation was plotted using 
representative values of ^ = 0.08 sec and k = 0.5* Comparison 
with tho thccrotical curves shews the sloj e of the power 
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spectral density of hole four to be approximtely the same as 
that for the thecrrtical curve. 



to 



COr^LUSICNS 



Ihe results wore inconclusive because the neutron break 
frequency was not determinable and as a result it was not 

1 -k 

possible to evaluate the ratio exactly. 

The best agreement betiveen theory and experimental 

results was obtained at z = 4 ft = 122 cm. At this position, 

the slope of the attenuated power spectral density function 

was approximately the same as that predicted by theory. The 

1 -k 

value of the ratio was determined to be approximately 
0.4/0.08 = 5«0. The results indicate that this ratio changes 
as z increases in value to beyond two diffusion lengths where 
It remains fairly constant until masked by the white noise 
of the detection chamber. 

The evaluation of the transfer function of the sub- 



critical assembly by analysis of the power spectral density 
shows promise and should be pursued further. 
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The plutorium-berylliim ncutrcn sources utilised in this 
investigation wore mounted on a pedest?.! beneath the assembly 
and were uximoderatod except for the moderation by the graphite 
of the assembly. Interesting results '^ight be obtained by 
encasing the sources in a block of paraffin moderator, and 
surrounding the moderated sources with graphite instead of 
the water shield, bore significant results night bo obtained 
by the use of a mechanical or electronic device for squaring 
and averaging tho data. This means of data processing would 
allow much inoro rapid and accurate detormln'^tion of the pov/er 
spectral density, hence broadening the scope of investigation 
possible. 
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Sample calrulaticns sr.ox/inr method 
spectral density. 

Hole three (z = 3 ft) 
/if = O.h cps, f = 0.6 cps 
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y = = 0,5669 









100 
= 0.323647 

^ ^ 0.638859 

= (yi~ y)^ = y^ - y^ . 

= 0.638859 - 0.3236U7 
= 0.315212 
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